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We present results from a survey carried out by the Balloon-borne Large Aper- 
ture Submillimeter Telescope (BLAST) on a 9 deg^ field near the South Ecliptic 
Pole at 250, 350 and 500 /im. The median la depths of the maps are 36.0, 26.4 
and 18.4 mJy, respectively. We apply a statistical method to estimate submil- 
limeter galaxy number counts and find that they are in agreement with other 
measurements made with the same instrument and with the more recent results 
from Herschel/SPIKE. Thanks to the large field observed, the new measurements 
give additional constraints on the bright end of the counts. We identify 132, 89 
and 61 sources with S/N > 4 at 250, 350, 500 /im, respectively and provide a 
multi-wavelength combined catalog of 232 sources with a significance > 4cr in 
at least one BLAST band. The new BLAST maps and catalogs are available 



publicly at http://blastexperiment.info 



Subject headings: cosmology: observations, submillimeter, galaxies: statistics, 
methods: data analysis 



1. Introduction 



Understanding the formation and evolution of galaxies is one of the foremost goals of 
experimental cosmology. In the redshift range z ~ 1—3, massive galaxies go through an evo- 
lutionary stage characterized by high rates of star formation. These early, dusty galaxies are 
best characterized by their thermal dust emission at far-infrared and submillimeter wave- 
length s. They are known to be the main comp onent of the Cosmic Infrared Background 
(CIB; iFixsen et all Il998l : iHauser fc Dwek 1 120011 ) . but still relatively httle is known about 
their nature and evolution. 

Observations b y the Balloon-borne Large Aperture Submillimeter Telescope (BLAST; 
Pascale et al.l 120081 ) have provided the first confusion-limited submillimeter maps at 250, 
350 and 500 /im, with a beam size of 36", 42" and 60", respectively, covering areas larger 
than 1 deg^. BLAST carried out surveys of the extragalactic sky in two blank fields, one 
center ed on the southern field of the Great Observatories Origins Deep Survey (GOODS- 
South; iDevlin et al.ll2009[ ) and one close to the South Ecliptic Pole (SEP). Each survey covers 
about lOdeg^. 

BLAST is the first instrument to provide maps of the sky at wavelengths near the 
peak of the CIB with enough sensitivity, sky coverage and angular resolution to identify 
a large number of sources, determine the detailed shape of the source counts and show 
that most of the FIRB comes from submillimeter sources already identified in deep 24 fun 
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several studies (JMarsden et al. 



12009 



surveys flDevlin et al.ll2009l: iPatanchon et al.l 12009 : iMarsden et al.ll2009l ). There have been 



Pascale et al.ll2009l ) of the statistical properties in the 



BLAST bands of sou rces from catalogs defined using Spitzer 24 /xm observations (FIDEL, 
Dickinson et al.l 120071 ). These have shown that 24yum sources may well contribute all of 
the CIB, with about half of the emission coming fr om galaxies at z >, 1.2 . Other statistica l 
analyses of the BLAST GOODS-South (BGS) data JPatanchon et al.ir2009l : Iviero et al.ll2009l ) 
have provided measurements of the differential number counts and of clustering on scales 
larger than the BLAST beam. 

In this paper we present results from the BLAST survey close to the SEP (hereafter 
BSEP field). This field was cho sen because of very lo w emission from infrared cirrus, inferred 
from the IRAS 100 /im map of ISchlegel et al.l fll998l ). Although point source detection (i.e. 
at high angular frequency) is not affected by fluctuations at large angular scales (of the 
order of a degree) from cirrus noise, a low-cirrus region is needed for studies of the CIB at 
these wavelengths. The SEP field does not have the same richness of multi-wavelength data 
as GOODS-South, but the multi-wavelength coverage is improving , in particular because 



of the observations carried out toward that region by the AKARI (JMatsuhara et al.ll2006 



Malek et al.ll2010l : iMatsuura et al.ll2010l ). AzTEC on the Atacama Submillimeter Telescope 
Experiment (Hatsukade et al. in prep.), the South Pole Telescope, the Atacama Cosmology 
Telescope and the Australia Telescope Compact Array. The catalog presented in this paper 
can provide submillimeter spectral energy distributions (^EDs) for the AzTEC and AKARI 
sources. Moreover, the BSEP field is one of the HerMEgil fields (ADF-S), so that the BSEP 
catalog can be used for cross-checking of the HerMES ADF-S data. 

We have performed an analysis of the statistical properties of the maps to estimate 
number counts and we have extracted sources using the Point Spread Funct ion (PSF) to 



filter the maps, applying the same techniques used in p revious BLAST pape rs (JDevhn et al. 



20091 : IPatanchon et al.ll2009l ). A more recent study by IChapin et al.l (120101 ) proposes a dif- 



ferent approach for filtering the submillimeter maps, namely a 'matched filter' optimized to 
maximize the signal-to-noise ratio (S/N) of individual point sources in the presence of mul- 
tiple noise sources. This approach significantly reduces the confusion noise, at the expense 
of slightly higher instrumental white noise. We have also used this approach for the BSEP 
data and have compared the results using both the PSF and the 'matched filter'. 



The same field has been also observed by Spitzer-MW^ at 24 and 70 /xm (IScott et al. 



20101). The combination of Spitzer and BLAST data enables the identification of mid-IR 
counterparts for ~ 50% of the BLAST sources and also provides useful limits for unidentified 



^http://hcrnics. Sussex. ac.uk 



-4- 



sources. The combined results can be used to constrain the SEDs of the galaxies across the 
entire Rayleigh- Jeans to Wien region of the far-IR spectrum and will be presented in a future 
paper. 

The paper is organized as follows: § 2 introduces the BLAST data; § 3 presents the 
statistical analysis used to calculate differential number counts and a comparison between 
number counts calculated in this field with previous measurements, as well as comparisons 
with the most recent evolutionary models; § 4 provides source lists constructed independently 
for each BLAST band and combined together and § 5 gives some brief conclusions. The 
source catalogs are provided in the Appendix. This paper constitutes the public release of 
the BLAST maps and multiwavelength catalogs of the BSEP field. 



2. BLAST Data 

During the BLAST 2006 Antarctic flight, deep observations (68 hours) were carried 
out over a 9deg field centered on (70.94°,— 53.50°) near the SEP. The BSEP maps are 
made from a large number of cross-linked scans, producing a uniform map with median la 
sensitivities (equivalent for detecting point sources) of 36.0mJy, 26.4 mJy, and 18.4 mJy, for 
the 250, 350, and 500 /im bands, respectively. Fig. [T]shows the combined signal-to-noise ratio 
image of all t hree BLAST bands . These new maps cover an area compa r able t o the BGS- 
Wide map of iDevlin et al.l (120091 ). We refer the reader to iPascale et al.l (120081) for a more 
detailed description of the characteristics of BLAST and iTruch et al.l ( l2009l ) for information 
on calibration and data reduction. 

The BLAST time-stream data were reduc ed using a cornmon pipeline to identify spikes, 
correct detector time drift and calibrate data (IPascale et al.ll2008t iTruch et al.ll2009l ). Maps 
were generated using the SANEPIC software, which uses a maximum-likelihood algorithm 
to estimate the optima l solution for the map, as well as producing an associated noise map 
( iPatanchon et al.ll2008l ). Absolute calibration is based on observations of the evolved star VY 



CMa and is estimat ed to have an unce rtainty of 10% (although strongly correlated between 
the three bands; see lTruch et al.l l2009l for details). 



While the maps represent the optimal weighting of the data across all spatial scales, 
the largest scales are less constrained due to various systematic effects. This can produce 
residual large-scale fluctuating patterns across the map. To suppress these spurious signals, 
all maps have been filtered to remove large scale angular frequencies, without affecting the 
contribution of individual point sources. This corresponds to scales in excess of about 10 
arcmin (approximately the size of the detector array projected on the sky). This procedure. 
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Fig. 1. — Map showing the combined signal-to- noise ratio of all three BLAST bands in the 
BSEP field. The circles mark the locations of the 232 sources detected with S/N > 4 in at 
least one band, as listed in Table IA4I The contours show the overlapping coverage of the 
11.5 deg^ Spitzer survey in this field (dark gray) and the region mapped at 90 fira by AKARI 
{light gray). 
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already used for the BGS maps (e.g. iDevlin et al.ll2009[ ). also explicitly sets the mean of each 



map to zero. The filtering induces negative shadows around the locations of bright sources in 
the maps, but this effect is fully taken into account in the analyses presented in this paper. 



The confusion noise in the map is estimated with the method explained in lMarsden et al. 



(120091 ). We fit a Gaussian to the distribution of the ratio of map pixel values to the instru- 
mental noise cXi. In the case of no confusion noise, the standard deviation of this Gaussian 
is unity by definition; if confusion noise is present, it gives the total noise of the map, at, in 
units of (Ti. The confusion noise, o"c, is then derived by subtracting the instrumental noise 
from the total noise in quadrature, (Xc = \/at — af. Results are reported in Table |5l As for 
the BGS- Wide map, the BSEP map is dominated by instrumental noise and so, using this 
criterion, the maps are not confusion limited. 

The BLAST maps in this release include three primary data products. The first set 
contains the 'raw' maps produced by SANEPIC. The second set of maps have been filtered 
spatially using a Wiener filter to whiten structures on scales larger than the instantaneous 
field of view. We also provide 'matched-filtered' maps (see § 14. 2p : these maps have been 
filtere d using a 'matched fi lter', which gives superior performs as maps approach the confusion 



limit fIChapin et all 120101 ) 



3. Number Counts 
3.1. P(D) analysis 



Number counts are estimated using the same method adopted in iPatanchon et al. 



(120091 ). In that paper it was shown that, in the S/N regime of BLAST, a statistical analy- 
sis of the maps is a better approach for obtaining number counts than counting individual 
sources. This is because it naturally allows for the correction of strong biases due to con- 
fusion and fiux boosting. This technique also has the advantage of providing an unbiased 
estimate of the counts at fiux densities well below the hmit at which sources can be detected 
individually. The method has been optimized to deal with inhomogeneous noise across the 
map and filtering to suppress large-scale noise. 

We start with a measure of the 'probability of defiection', P(D), the histogram of 
pi xel values, and use thi s to try to obtain the best estimate of the number counts. As 



m 



Patanchon et al.l ( 12009! ) , we choose to parametrize the differential number counts by a set 
of amplitudes at a few predefined fiuxes, with the intervals between fiux nodes interpolated 
with power-laws to impose continuity of the counts. The number and locations of the nodes 
are chosen as a compromise between increasing the number of free parameters, to give a 
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better representation of the true underlying counts, and keeping the number of nodes small 
enough to provide useful estimates of their values. We find that no more than about four 
amplitude parameters can be estimated for each waveband. This is lower than the numbers 
estimated from the BGS maps (up to seven nodes at 250 /xm), but this is not unexpected 
given the extra dynamic range of BGS- Wide plus BGS-Deep. Our observations provide useful 
information starting from ~ 0.05 Jy and brighter. A larger number of nodes would increase 
the correlation between neighbouring nodes and thus would not add any information. 

The P(D) analysis is carried out by minimizing the negative log-likelihood 

^9) = -J2n^log{pi{9)) - log(iV!) + ^log(n,!), (1) 

i i 

where rii is the number of pixels with flux densities in the ith flux bin interval; Pi is the result 
of the integral of the P(D) in the ith bin, normalized such that YlPi — Ij N is the total 
number of measurements (pixels in the map); 6 are the parameters of the number counts 
model and the last two terms derive from the normalization of the multinomial distribution 



function (see iPatanchon et al.l l2009l for a derivation). Eq. [T] is only strictly correct if the 
probability distribution is the same for all observations (pixels). This is not the case for 
the BGS maps, with the expected distributions being different for the deep and the wide 
regions, and it is not completely true for the BSEP maps, even though differences in the 
noise levels are not as dramatic as in the BGS case. To deal with this we have divided the 
observed maps into a small number of regions (4), such that in each zone the noise variance 
is approximately constant. We then compute the quantities in Eq. [T]for each of the regions. 

The P(D) analysis works under the hypothesis that all the sources are point sources, 
i.e. all sources in the map have the size of the beam. This is not exactly true in our BSEP 
maps, where at least one local galaxy is partly resolved by the instrument and appears as an 
extended source (see § 14. 3p . Before performing the fit, we have masked the extended source 
and its dark ring. 

We have chosen to fit power-laws for differential number counts within predefined flux 
density bins, as described above. Five distinct power-laws are estimated (a total of six 
free parameters) for the differential number counts at 250 ^m and four power-laws (five 
parameters) at both 350 and 500 /xm. The choice of flux densities for the boundary nodes 
is set by requiring them to be very far from the typical values constrained by BLAST, 
so that the result is independent of our particular choice. Best-fit number counts for the 
three wavelengths are presented in Fig. |2] and Table [U The uncertainties shown in the 
figure are computed from the 68% confidence intervals on the marginalized distributions for 
each parameter separately. The marginal distributions (see Fig. |3]) have been estimated by 
sampling the likelihood around its maximum using a Markov Chain Monte Carlo Metropolis 



Hastings method (MCMCMH; IChib fc Greenberg Ill995l ). The median values presented in 
Table [Dare not exactly the parameters of the best-fit model, due to non-Gaussian likelihoods 
around the maximum. Pearson correlation matrices for the parameters are given in Tables [H 
ElandU 

We compare the predicted histograms of the best-fit multi-power-law model with the 
actual histograms of the maps in Fig. HI We plot the histograms of the two zones that include 
about ~ 95% of the pixels and ignore the others here, since they only give weak additional 
constraints on the parameters. 



Patanchon et al.l (J2009[ ) have demonstrated that, in a noise-dominated regime, it is useful 



to cross-correlate the maps with the beam kernel before P(D) analysis, even if this might 
not be the optimal choice. We have performed the analysis both convolving the maps with 
the same instrumental beam already used for the BGS maps and using a 'matched filter', 
optimized for confusion-hmited maps (see § 14. 2p . As already pointed out in § |2l the BSEP 
maps are dominated by instrumental noise, thus the 'matched filter' does not significantly 
improve the analysis: the differences between the two methods are much smaller that the 
uncertainties on the counts. 

The technique used in this paper also assumes that galaxies are randomly and inde- 
pendently distributed over the sky. However, sign ificant correlations have been found in the 
background of the BGS maps f IViero et al.l 120091 ), mainly due to cl u sterin g on scales larger 
than the BLAST beam. As already discussed in iPatanchon et al.l (12009! ) . the influence of 
clustering on the measurements of number counts is negligible with respect to the uncertain- 
ties. The clustering in BSEP is comparable to what is m easured in BG S, despite having the 
outskirts of a local cluster (DRCG 0428—53 at z = 0.04, iDressel Ill980l ) cover almost half of 
the map. As a test, we repeated the analysis on a sub-map that is free of cluster members. 
Our previous study on a galaxy cluster suggests that the contribution of the cluster to the 
submill imeter number coun ts may be significant in particular at the faint-end flux node in 
BSEP f lBraglia et al.ll2010l ). but the counts derived from the sub- map are consistent with 
those derived from the full map at the faint end of the distribution. At the bright end the 
counts from the sub-map are lower, but consistent within the small number statistics. In 
conclusion, the influence of the cluster can be considered negligible. 



3.2. Comparison with other measurements 

The number counts provided in this work can be compared with previously published 
counts in the same bands. All the existing differential number counts come from very recent 
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Fig. 2. — Best-fit differential number counts at the three BLAST wavelengths {solid lines). 
Uncertainties are derived from the marginalized distribution for each parameter {dotted 
lines); because of the non-Gaussian behavior of the likelihood around its maximum, the 
best-fit model is not centered on the error bars. The first and last power- laws {dotted lines) 
are upper limits . Number counts in the same bands fr om other st udies are also p l otted : 



Patanchon et all (120091 ) {diamonds); lOliver et all (120 lOl ) {squares); I Clements et al 

{crosses). 
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Table 1. Best-fit differential number counts 





250 Aim 






350 ^m 






500 fim 




Node 


Best fit 


Marginal 


Node 


Best fit 


Marginal 


Node 


Best fit 


Marginal 


[Jy] 


[log(deg 


:^Jy-^)] 


[Jy] 


[log(de 


g'Jy-')] 


[Jy] 


[log(de 


g^Jy-^)] 


10-4 


9.39 


< 10.09 


5 X 10-5 


11.64 


< 12.10 


2.5 X 10-5 


12.43 


< 12.81 


0.1 


3.23 


q 97+0.11 
^■^'-0.10 


0.05 


3.86 


3.89lg;l? 


0.03 


4.08 


4 24+0-13 


0.25 


1.01 


OMt'oit 


0.15 


1.16 


0.93l°:i 


0.08 


1.66 


1.24t?l4 


0.5 


0.65 


0-751:^21 


0.5 


0.31 


n ^1+0-23 
^■"^-0.28 


0.25 


0.58 


r. f.7+0.32 

•^■"'-0.38 


1.2 


-0.080 


p. 9r.+0.22 


5 


-4.09 


< -2.24 


2.5 


-11.17 


< -4.27 


10 


-14.42 


< -6.78 















Table 2. Pearson correlation matrix for the parameterized dN/dS model at 250 /im 
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Note. — Coefficients are computed as the covariance of the 
distributions of the parameters around the maximum of the likeli- 
hood. The distributions were derived using MCMCMH method. 
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Table 3. Pearson correlation matrix for the parameterized dN/dS model at 350 /xm 
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Note. — Coefficients are computed as the covariance of 
tlie distributions of the parameters around the maximum of 
the hkehhood. The distributions were derived using MCM- 
CMH method. 



Table 4. Pearson correlation matrix for the parameterized dN/dS model at 500 /xm 
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Note. — Coefficients are computed as the covariance of the 
distributions of the parameters around the maximum of the 
hkehhood. The distributions were derived using MCMCMH 
method. 
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Fig. 3. — Likelihood distributions and contours for pairs of parameters associated with 
adjacent nodes at 250 /xm. The two curves in each panel represent 68% and 95% intervals 
and the solid and dashed lines represent the median and la dispersion for each marginalized 
parameter: the distributions of the parameters are estimated from sampling the likelihood 
with MCMCMH. The red crosses mark the best-fit values. 
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Fig. 4. — Histograms of pixel values for two regions having different noise level, including 
about ~ 95% of pixels of the BSEP map, compared with predictions of the best-fit model 
of the differential counts at 250 /im. The corresponding plots at 350 and 500 fira look very 
similar. 
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maps provided by BLAST and Herschel/SPIKE. In Fig. [2] we com pare our best-fit e s timat e 
for the diffe r ential nur nber counts with the published c ounts by iPatanchon et al.l ( 120091 ) 



Oliver et all (120101 ) and lClements et all (120101) . While the 



Patanchon et al. 



derived from a statistical analysis of the map, lOliver et al. 



(l2009h counts are 



mm and lClements et al.l tom 



extract candidate sources and apply corrections to the raw number counts for Eddington 
bias, reliability and incompleteness. Despite these differences, the number counts from these 
studies are in very good agreement with our results. 

The sky coverage delivered by the BSEP map provides further constraints on the number 
counts at the bright end of the distribution. 



3.3. Comparison with models 



Fi g. \5\ compares our be st-fit counts with the model of IValiante et al.l (120091 ). As pointed 
out by lOliver et al.l (120101 ). this model has the distinct advantage, compared with most 
published models, of being able to fit the rise in the counts for fiuxes in the range 10—100 mJy. 
This could suggest that galaxies at high redshift have on average a colder temperature than 
local objects with the same luminosity. However, with suitable changes in the luminosity 
functions, it may be possible to make other models fit the observed counts, so it is premature 
to conclude that the spectral energy distributions need revision. Additionally, the model is 
able to reproduce the measured counts within the error bars in the range 50 — 1000 mJy, as 
explored in these new BSEP data. 



4. Source Catalogs 

We have compiled a catalog of point sources with fiux S > 3a for each band (see Ta- 
bles lAll IA2[ IA3p using a source-finding algorithm which selects the peaks in a smoothed 
map produced by the noise-weighted convolution of the image with the telescope PSF 
(JDevlin et al.ll2009l ). The adopted a is the total noise of the map, o"t, defined in § |2l The 
source lists were synthesized into a common catalog using a procedure which accounts for 
the significance and positional uncertainty of the counterparts in each band. The radius of 
the la positional error circle, cxp, for a submillimeter galaxy in a catalog with signal to noise 
/i which has not been corrected for the Eddington bias type of fiux boosting can be expressed 
as 



ap = 0.99[fi^ - {2a + 4)] 



-1/2 



(2) 
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Fig. 5. — Best-fit differential number count s for the three BLAS T wavelengths (circles) 
compared with a reahzation of the model of IValiante et al.l (120091 ) [solid line). The plots 
also show the cruder estimate of the number counts derived directly from the > 4a catalogs 
after attempting to correct for completeness and false detection rate (triangles) . The vertical 
lines define the flux limits of the catalogs. Note that the error bars are not marginalized in 
the same way as for the P(D) counts, so that it is not possible to make a simple comparison. 
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for power-law counts of the form N(> S) oc S °, wh ere 6 is the full width at half maximum 



(FWHM) of the telescope beam ( llvison et al.l 120071 ). Using this formula, error circles were 



calculated assuming; a = 2, as obtained by fitting the number counts with a single power-law 



( iPatanchon et al.ll2009[ ). A minimum la error circle of 5" was imposed, equal to the la 



pointing uncertainty of the maps. 

The combined catalog is comprised of all sources with a significance > Aa in at least 
one band. Sources from other bands are considered to be matches if they are located within 
twice the radius of their respective error circles added in quadrature. Positions of sources in 
the resulting combined catalog were computed by averaging all the positions, weighted by 

The number of sources detected at > Aa by band is 132, 89 and 61 at 250, 350 and 
500 /im, respectively (see Table E]). Table R4] lists the coordinates and 250, 350 and 500 /xm 
flux densities and uncertainties for the 232 sources of the combined catalog. 



4.1. Completeness and false detection rates 

'Completeness' is the probability of detecting a source of a given intrinsic flux density, 
given the depth of the survey and the source extraction algorithm. We estimate completeness 
by adding to the initial maps 1000 artificial point sources at random positions and calculating 
their recovery rates after performing the same source extraction procedure as for the real 
data. The artificial point sources are modeled as the PSF scaled by the flux density. To 
avoid blending the simulated sources together or with the existing sources, every simulated 
source is injected at a distance larger than 2.5 times the half width at half maximum from 
any real sources and from any other simulated ones. In principle this procedure introduces 
a bias which we estimate to be well under 1%. The number density of artificial sources is 
such that it does not appreciably change the noise properties of the maps. A source was 
considered to be detected if there was a detection within a circle centered on the source and 
within the radius of the FWHM of the corresponding wavelength. The catalog completeness 
as a function of intrinsic flux density is shown in Fig. [61 Table |5] gives the 50%, 80% and 
95% completeness flux densities at each wavelength. 

We calculate the false detection rates (FDRs) by performing the extraction algorithm 
on 100 simulated noise realization maps. This is a c onservative oyerest imate of the FDRs 



and has the advantage of being model-independent (jPerera et al.ll2008l ). FDR values are 
reported in Table [5] for sources with a significance of > 4a. At fluxes > 5a they are all 
consistent with zero. 
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Fig. [5] shows the number counts derived from the > 4a catalogs after correcting for 
completeness and false detection rate. We can see that the catalog provides counts in good 
reasonably agreement with those derived by the P(D) analysis, but for a restricted flux 
density range. 



4.2. Matched filter 

We have repeated the source extraction procedure along with completeness and false 
detection rate tests on ma ps which have been c onvolved with the 'matched filter', calculated 



following the algorithm of IChapin et al.l (120 lOf ) instead of the PSF. This filter is optimized 



to maximize the S/N of individual point sources in the presence of noise: it significantly 
reduces the confusion noise and slightly increases the white instrumental noise. This filter 
should perform significantly better as maps become more dominated by confusion. 

Since the BSEP maps are dominated by instrumental noise, we find little benefit in the 
use of the 'matched filter', as shown in Table |6l Even though we identify a larger number of 
detected sources using the 'matched filter', the fraction of false detections is also larger. For 
this reason, we prefer to publish the more conservative, and easily reproducible, source lists 
derived by filtering the maps with the PSF. 



4.3. Extended sources 

The technique we have used to estimate the fiux densities of the BLAST sources is not 
accurate in the case of extended sources. In the BSEP map there is one clearly extended 
source, identified with the galaxy NGC 1617 and detected as BLAST J043137-543604 in 
our catalog. We have performed aperture photometry on this source and have estimated the 
uncertainty on the fiux density taking into account both instrumental noise and calibration 
uncertainties. The measured values are reported in Table IA4[ 



5. Conclusions 

We used BLAST to image a 9 deg^ field near the SEP at 250, 350 and 500 /im, achieving 
median la depths of 36.0, 26.4 and 18.4 mJy at each wavelength. We have identified 132, 
89 and 61 sources with S/N > Aa at 250, 350 and 500 /zm, respectively. These have been 
compiled into a multi-wavelength catalog of 232 sources with a significance > Aa in at least 
one BLAST band. 
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Using the P(D) technique, we have measured the differential number counts up to 1.2, 
0.5 and 0.25 Jy at 250, 350 and 500 /xm, respectively. The new measurements agree with 
previous results from BLAST and more recent results from SPIRE and, thanks to the large 
area observed, give improved constrains at the bright end of the counts. 



We have released the BLAST maps and catalogs to the public at http://blastexperiment.info 
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Table 5. BLAST-SEP maps and catalogs 



Band ai ctc Number False 50% 80% yo7o 

of sources detections completeness completeness completeness 

[yum] [mJy] [mJy] > 4a > Aa [mJy] [mJy] [mJy] 



250 


36.0 


20.6 


132 


2.28 


192 


241 


299 


350 


26.4 


18.2 


89 


0.94 


137 


179 


220 


500 


18.4 


15.2 


61 


0.14 


97 


132 


157 



-19- 



Table 6. BLAST-SEP maps and catalogs with Matched Filter 



Band a-, a^ Number False 50% 80% 95% 

of sources detections completeness completeness completeness 

[;um] [mJy] [mJy] > Aa > Act [mJy] [mJy] [mJy] 



250 


38.2 


16.5 


145 


7.53 


194 


241 


297 


350 


28.9 


14.2 


100 


6.80 


144 


184 


226 


500 


20.7 


11.6 


90 


3.13 


99 


134 


163 



E 0.4 




t 0,4 





E 0.4 - 



.2 0.3 0.4 

Flux [Jy] 



Fig. 6. — Catalog completeness at 250, 350 and 500 /xm. The error bars are estimated from 
Icr binomial uncertainties. 
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7. Appendix A: data tables 

The following data tables are provided in the Appendix: the lists of subniillimeter 
sources with flux > 3a at 250 /im (Table |AT]), 350 /xm (Table |M]) and 500 /im (Table IMJ; 
the multi-wavelength combined catalog of sources with a signiflcance > 4cr in at least one 
band (Table R4|) . Data tables are published in their entirety in the electronic edition of the 
Astrophysical Journal Supplement Series. The flrst 25 entries are shown here for guidance 
regarding their form and content. 
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Table Al. Catalog of BLAST 250 yum sources in SEP with significance > 3a. Flux 
densities are not corrected for Eddington-type bias. 



ID 



BLAST ID 



RA 

(J2000) 



Dec 

(J2000) 



'5'250 

Jy 



17250 

Jy 



S/N 



1 


BLAST250 J043138-543604 


67.909134 


-54.601269 


2.O51 


O.211 


9.76 


2 


BLAST250 J043516-541902 


68.818535 


-54.317352 


0.818 


0.041 


19.77 


3 


BLAST250 J043126-542507 


67.860474 


-54.418732 


0.698 


0.043 


16.15 


4 


BLAST250 J043422-535358 


68.591782 


-53.899590 


0.549 


0.041 


13.40 


5 


BLAST250 J044936-535427 


72.400566 


-53.907761 


0.515 


0.040 


12.78 


6 


BLAST250 J042755-550332 


66.980072 


-55.058910 


0.716 


0.053 


13.61 


7 


BLAST250 J045412-532127 


73.550888 


-53.357536 


0.482 


0.041 


11.78 


8 


BLAST250 J045053-531233 


72.724892 


-53.209232 


0.482 


0.042 


11.56 


9 


BLAST250 J043424-544132 


68.601288 


-54.692341 


0.478 


0.042 


11.44 


10 


BLAST250 J044801-532609 


72.005684 


-53.435989 


0.467 


0.043 


10.97 


11 


BLAST250 J045045-533144 


72.687759 


-53.528919 


0.416 


0.040 


10.42 


12 


BLAST250 J045624-523057 


74.101578 


-52.515938 


0.361 


0.040 


9.10 


13 


BLAST250 J045313-524746 


73.306427 


-52.796268 


0.345 


0.040 


8.69 


14 


BLAST250 J044839-535429 


72.165001 


-53.908199 


0.339 


0.041 


8.36 


15 


BLAST250 J045307-525416 


73.282074 


-52.904690 


0.337 


0.040 


8.34 


16 


BLAST250 J042705-541554 


66.773399 


-54.265018 


0.333 


0.040 


8.23 


17 


BLAST250 J044706-524028 


71.777534 


-52.674500 


0.320 


0.039 


8.10 


18 


BLAST250 J045759-524129 


74.499512 


-52.691544 


0.331 


0.041 


8.18 


19 


BLAST250 J043553~542417 


68.973793 


-54.404816 


0.351 


0.042 


8.28 


20 


BLAST250 J043513-541230 


68.807793 


-54.208553 


0.329 


0.041 


8.04 


21 


BLAST250 J045443-530540 


73.680298 


-53.094635 


0.314 


0.041 


7.63 


22 


BLAST250 J045139-534459 


72.916176 


-53.749763 


0.305 


0.040 


7.55 


23 


BLAST250 J044853-523039 


72.221237 


-52.511059 


0.306 


0.041 


7.53 


24 


BLAST250 J044125-532834 


70.357887 


-53.476337 


0.318 


0.042 


7.59 


25 


BLAST250 J043346-534651 


68.441719 


-53.781109 


0.295 


0.040 


7.38 



These flux densities come from aperture photometry. See § 14.31 for details. 
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Table A2. Catalog of BLAST 350 yum sources in SEP with significance > 3a. Flux 
densities are not corrected for Eddington-type bias. 



ID 


BLAST ID 


RA 


Dec 


'S'350 


o'sso 


S/N 






(J2000) 


(J2000) 


Jy 


Jy 




1 


BLAST350 J043137-543603 


67.904358 


-54.601097 


1.091 


0.121 


9.08 


2 


BLAST350 J043516-541902 


68.818398 


-54.317253 


0.351 


0.032 


10.92 


3 


BLAST350 J042754-550331 


66.975388 


"55.058666 


0.472 


0.044 


10.80 


4 


BLAST350 J043125-542516 


67.855537 


-54.421322 


0.325 


0.033 


9.88 


5 


BLAST350 J043421-535408 


68.591515 


-53.902260 


0.280 


0.031 


8.93 


6 


BLAST350 J045412-532127 


73.550613 


-53.357510 


0.272 


0.031 


8.66 


7 


BLAST350 J044935-535437 


72.396095 


-53.910458 


0.248 


0.032 


7.80 


8 


BLAST350 J045044-533143 


72.687477 


-53.528831 


0.229 


0.031 


7.37 


9 


BLAST350 J044706-524018 


71.777740 


-52.671692 


0.204 


0.031 


6.60 


10 


BLAST350 J043424-544132 


68.601440 


-54.692429 


0.198 


0.032 


6.11 


11 


BLAST350 J045317-533816 


73.321846 


-53.637840 


0.185 


0.031 


5.87 


12 


BLAST350 J045053-531233 


72.724693 


-53.209290 


0.197 


0.033 


5.97 


13 


BLAST350 J045025~524125 


72.607147 


-52.690285 


0.180 


0.031 


5.77 


14 


BLAST350 J043624-542508 


69.101585 


-54.419037 


0.178 


0.033 


5.32 


15 


BLAST350 J043935-530926 


69.898155 


-53.157463 


0.217 


0.039 


5.53 


16 


BLAST350 J044839-535428 


72.164772 


-53.908028 


0.152 


0.030 


5.00 


17 


BLAST350 J044436-540946 


71.150352 


-54.162781 


0.165 


0.032 


5.10 


18 


BLAST350 J044853-523039 


72.221313 


-52.510956 


0.157 


0.031 


5.01 


19 


BLAST350 J044802-532608 


72.010269 


-53.435829 


0.172 


0.033 


5.14 


20 


BLAST350 J044040-542055 


70.169220 


-54.348644 


0.161 


0.032 


4.96 


21 


BLAST350 J045623-523106 


74.097374 


-52.518547 


0.148 


0.031 


4.80 


22 


BLAST350 J044944-525427 


72.434944 


-52.907745 


0.146 


0.031 


4.77 


23 


BLAST350 J044831-540229 


72.132347 


-54.041489 


0.154 


0.032 


4.85 


24 


BLAST350 J043512-541210 


68.803497 


-54.202976 


0.153 


0.032 


4.81 


25 


BLAST350 J043830-541840 


69.628357 


-54.311268 


0.153 


0.032 


4.77 



These flux densities come from aperture photometry. See § 14.31 for details. 
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Table A3. Catalog of BLAST 500 yum sources in SEP with significance > 3a. Flux 
densities are not corrected for Eddington-type bias. 



ID 



BLAST ID 



RA 

(J2000) 



Dec 

(J2000) 



S500 

Jy 



fsoo 

Jy 



S/N 



1 


BLAST500 J043137-543603 


67.904343 


-54.601093 


O.3951 


O.O551 


7.18 


2 


BLAST500 J042756-550322 


66.985207 


-55.056328 


0.257 


0.030 


8.52 


3 


BLAST500 J043125-542506 


67.855965 


-54.418583 


0.167 


0.025 


6.73 


4 


BLAST500 J045056-531633 


72.735146 


-53.275887 


0.150 


0.023 


6.40 


5 


BLAST500 J043516-541912 


68.818130 


-54.320080 


0.147 


0.024 


6.21 


6 


BLAST500 J043521-550056 


68.840088 


-55.015621 


0.147 


0.026 


5.73 


7 


BLAST500 J044742-533019 


71.926003 


-53.505398 


0.130 


0.024 


5.36 


8 


BLAST500 J045600-524902 


74.001236 


-52.817482 


0.115 


0.023 


5.04 


9 


BLAST500 J042640-541057 


66.669762 


-54.182514 


0.150 


0.027 


5.56 


10 


BLAST500 J045046-533133 


72.692162 


-53.526104 


0.114 


0.023 


5.00 


11 


BLAST500 J045307-525357 


73.281929 


-52.899330 


0.118 


0.023 


5.06 


12 


BLAST500 J050244-525150 


75.686989 


-52.864151 


0.154 


0.028 


5.57 


13 


BLAST500 J045920-521402 


74.836861 


-52.234043 


0.116 


0.023 


4.96 


14 


BLAST500 J045107-525202 


72.779480 


-52.867264 


0.109 


0.023 


4.78 


15 


BLAST500 J044156-531041 


70.485947 


-53.178085 


0.115 


0.024 


4.82 


16 


BLAST500 J050215-523923 


75.564407 


-52.656528 


0.130 


0.026 


5.04 


17 


BLAST500 J043750-532126 


69.459808 


-53.357227 


0.120 


0.025 


4.89 


18 


BLAST500 J044350-525220 


70.961983 


-52.872311 


0.118 


0.024 


4.85 


19 


BLAST500 J043412-545017 


68.550507 


-54.838169 


0.121 


0.025 


4.89 


20 


BLAST500 J044310-542034 


70.793472 


-54.343006 


0.118 


0.024 


4.84 


21 


BLAST500 J045840-523546 


74.668068 


-52.596275 


0.120 


0.025 


4.83 


22 


BLAST500 J045442-530530 


73.675896 


-53.091942 


0.110 


0.024 


4.63 


23 


BLAST500 J044154-540351 


70.475601 


-54.064274 


0.104 


0.023 


4.49 


24 


BLAST500 J044232-535158 


70.633949 


-53.866184 


0.106 


0.024 


4.50 


25 


BLAST500 J043832-541811 


69.633400 


-54.303066 


0.107 


0.024 


4.48 



These flux densities come from aperture photometry. See § 14.31 for details. 
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Table A4. Multi-wavelength catalog of BLAST sources in SEP with > 4a detection in at 

least one band. Flux densities are not corrected for Eddington-type bias. The 

non-detections are expressed as 3a upper limits. 



ID 



BLAST ID 



RA 

(J2000) 



Dec 
(J2000) 



»S'250 

Jy 



f250 

Jy 



S350 
Jy 



Jy 



S5QQ 
Jy 



Jy 



1 


BLAST J043137-543604 


67.905945 


-54.601154 


2.O51 


O.211 


I.O91 


O.I21 


O.3951 


O.O551 


2 


BLAST J043516-541904 


68.818398 


-54.317852 


0.818 


0.041 


0.351 


0.032 


0.147 


0.024 


3 


BLAST J043125-542510 


67.857567 


-54.419716 


0.698 


0.043 


0.325 


0.033 


0.167 


0.025 


4 


BLAST J043422-535403 


68.591995 


-53.901043 


0.549 


0.041 


0.280 


0.031 


0.094 


0.023 


5 


BLAST J044935-535432 


72.398254 


-53.909065 


0.515 


0.040 


0.248 


0.032 


0.070 


0.023 


6 


BLAST J042755-550329 


66.979797 


-55.058105 


0.716 


0.053 


0.472 


0.044 


0.257 


0.030 


7 


BLAST J045412-532127 


73.550751 


-53.357513 


0.482 


0.041 


0.272 


0.031 


0.092 


0.024 


8 


BLAST J045053-531232 


72.724510 


-53.208916 


0.482 


0.042 


0.197 


0.033 


0.085 


0.024 


9 


BLAST J043424-544132 


68.601357 


-54.692379 


0.478 


0.042 


0.198 


0.032 


<0.074 


- 


10 


BLAST J044801-532609 


72.007538 


-53.435932 


0.467 


0.043 


0.172 


0.033 


0.097 


0.024 


11 


BLAST J045045-533142 


72.688240 


-53.528496 


0.416 


0.040 


0.229 


0.031 


0.114 


0.023 


12 


BLAST J045624~523100 


74.100220 


-52.516781 


0.361 


0.040 


0.148 


0.031 


<0.071 


- 


13 


BLAST J045313-524746 


73.306450 


-52.796272 


0.345 


0.040 


0.129 


0.031 


<0.069 


- 


14 


BLAST J044839-535429 


72.164917 


-53.908138 


0.339 


0.041 


0.152 


0.030 


<0.069 


- 


15 


BLAST J045307-525411 


73.282059 


-52.903137 


0.337 


0.040 


0.137 


0.031 


0.118 


0.023 


16 


BLAST J042705-541554 


66.773399 


-54.265018 


0.333 


0.040 


<0.095 


- 


<0.072 


- 


17 


BLAST J044706-524023 


71.777649 


-52.673302 


0.320 


0.039 


0.204 


0.031 


0.073 


0.022 


18 


BLAST J045759-524129 


74.499466 


-52.691570 


0.331 


0.041 


0.131 


0.032 


<0.069 


- 


19 


BLAST J043553-542417 


68.973129 


-54.404781 


0.351 


0.042 


0.117 


0.033 


<0.072 


- 


20 


BLAST J043513-541228 


68.807434 


-54.207802 


0.329 


0.041 


<0.095 


- 


0.081 


0.024 


21 


BLAST J045443~530537 


73.679634 


-53.093685 


0.314 


0.041 


0.132 


0.032 


0.110 


0.024 


22 


BLAST J045139-534501 


72.916527 


-53.750340 


0.305 


0.040 


<0.094 


- 


0.073 


0.023 


23 


BLAST J044853-523039 


72.221260 


-52.511024 


0.306 


0.041 


0.157 


0.031 


<0.071 


- 


24 


BLAST J044126-532832 


70.359184 


-53.475574 


0.318 


0.042 


0.146 


0.032 


<0.073 


- 


25 


BLAST J043346-534651 


68.441719 


-53.781109 


0.295 


0.040 


<0.096 


- 


<0.073 


- 



These flux densities come from aperture photometry. See ij 14.31 for details. 



